Decoherence strength of multiple non-Markovian environments 
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It is known that one can characterize the decoherence strength of a Markovian environment by the 
product of its temperature and induced damping, and order the decoherence strength of multiple 
environments by this quantity. We show that for non-Markovian environments in the weak coupling 
regime there also exists a natural (albeit partial) ordering of environment-induced irreversibility 
within a perturbative treatment. This measure can be applied to both low-temperature and non- 
equilibrium environments. 



I. INTRODUCTION 

Environment-induced decoherence is an essential pro- 
cess for a quantum system to acquire classical attributes 
[1-3]. To characterize how strong an environment can 
induce decoherence in an open quantum system it is 
desirable to come up with a measure of the "decoher- 
ence strength" of each environment acting on the sys- 
tem. For a Markovian environment, which for quantum 
Brownian motion refers to high temperatures and ohmic 
spectral density, such a measure can be constructed by 
the product of its temperature and damping rate. How- 
ever, such a measure may not exist for a general environ- 
ment (with non-ohmic spectral density functions and un- 
der low-temperature conditions; see or, even more 
challenging, for nonequilibrium environments, where the 
notion of temperature loses meaning. 

We show in this paper that at least perturbatively 
for weak coupling between the system and these envi- 
ronments, low-temperature and non-equilibrium environ- 
ments can be partially ordered. The object of compar- 
ison is, in fact, the correlation function of the collective 
coupling operators of the environment. As with matrices 
and kernels, quantum correlations can only be partially 
ordered. However, this does not rule out nontrivial com- 
parisons: in Sec. Ill we explicitly detail how one can com- 
pare decoherence strengths for combinations of thermal 
reservoirs without resorting to the concocted notion of 
an effective temperature, since in general it does not ex- 
ist. Our general relation includes some recently reported 
results [7] as special cases. 

It should be noted that strictly speaking it is the order- 
ing of environment-induced irreversibility (in the sense 
of contraction of the system's state-space volume un- 
der time evolution) which will be established here. This 
involves the phenomena of decoherence, typically domi- 
nated by information flow to the environment, as well as 
thermalization (or equilibration), where energy flow also 
plays a crucial role. There are many interesting situa- 
tions (especially when the system-environment coupling 
is not very strong) in which both phenomena correspond 




FIG. 1. Consider, for instance, a system of nano- mechanical 
resonantors Sm interacting with a system of optical modes So 
3—10] , wherein the optical modes experience dissipation and 
thermal noise To from the cavity field yet the resonators ex- 
perience dissipation and thermal noise Tm from a phonon en- 
vironment. The combined environment is out of equilibrium 
and cannot be described by a single spectral-density func- 
tion and temperature since it does not obey the fluctuation- 
dissipation relation. Therefore, the decoherence of this multi- 
partite system cannot be characterized in terms of a temper- 
ature and a dissipation coefficient. 



to rather different timescales and one can focus on deco- 
herence as the source of irreversibility at short times. 
Throughout the paper we use units with h = 1. 



II. QUANTUM CORRELATIONS & 
DECOHERENCE STRENGTH 

Consider a quantum system weakly interacting with 
an environment and with the interaction Hamiltonian ex- 
panded as a sum of separable operators: 

Hi = ^ L„ (g) 1„ , (1) 

n 
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where L„ and 1„ are system and environment opera- 
tors respectively. [The interaction Hamiltonian as well 
as the set of operators L n and \ n in Eq. (1) are in the 
Schrodinger picture, but they can in general be time de- 
pendent.] 

The environment coupling operators \ n will typically 
be collective observables of the environment, with depen- 
dence upon very many modes. The system-environment 
interaction will be treated perturbatively, so that the cen- 
tral ingredient is the (multivariate) correlation function 
of the environment: 

a nm (t,T) = {l n (t)l m (T)) E , (2) 

where l n (t) represents the time-evolving 1„ in the inter- 
action (Dirac) picture. In the influence functional for- 
malism [11] for the quantum Brownian model with bi- 
linear couplings between the system and its environment 
[4, 5, 12] the correlation function appears as the kernel 
in the exponent of a Gaussian influence functional, called 
the influence kernel £ in Refs. [13, 14]. Alternatively, in 
quantum state diffusion [15] this kernel takes the explicit 
role of a correlator for complex Gaussian noise. The in- 
fluence kernel, or equivalently, the complex correlation 
function, can be written as a sum of two real parts cor- 
responding to the noise and dissipation kernels [13, 14]: 

a(t,r) =v(Vr)+t m(*,t) . (3) 

complex noise noise dissipation 

The noise kernel v appears within the influence func- 
tional formalism as the correlation of an ordinary real 
stochastic source, whereas the dissipation kernel p alone 
would produce a purely homogeneous (though not neces- 
sarily positivity preserving) evolution. These same roles 
can also be inferred from the Heisenberg equations of 
motion for the system operators after integrating the en- 
vironment dynamics, producing the so-called quantum 
Langevin equation [16]. 

For weak coupling one can treat perturbatively the 
system-environment interaction. Using the notation of 
Ref. [ ], the second-order master equation [18-20] for 
the reduced density matrix p of the system can be rep- 
resented in terms of the noise correlation as 

p=[-iU,p]+C 2 {p}, (4) 

with the second-order contribution given by 

£-2{p\ = ^ [ L,i ' P ( A «™ <> L ™) f - (A„ m o L m ) p] , 

nm 

(5) 

where the A operators and o product define the second- 
order operators 



(A nm oL m )(<) = 




(t,T) {g (t,T)L m (T)} , (6) 



with L m (r) in the Schrodinger picture and the superop- 
erator Qo{t,r) being the free system propagator, which 



for a time-independent Hamiltonian simply corresponds 
to Qo{t,T)p(r) = e - l ( t ^ 7 ") H p(r)e l ( t - r ) H . Because of the 
nonlocal character of the noise correlation, the time- 
translation generator in Eq. (5) is not generally of Lind- 
blad form. The theorem by Lindblad [21] and Gorini, 
Kossakowski and Sudarshan (GKS) [22] specifically char- 
acterizes the algebraic generators $ for all completely- 
positive maps e v * (77 > parameterizes the semi-group 
[21, 22] in the Markovian case and will be taken equal to 
one here) as 

*P= -i[®,p]+^2A IJ (e I pe\- \{ e j e i>p}) ' 

unitary v j 
decohcrent 

( ? ) 

where ej = for some basis {]«)} of the system's 

Hilbert space and with I = ii' labeling all possible pairs 
of indices. Such generators and the dynamics they engen- 
der when the master equation has the Lindblad form have 
been extensively studied [23-32]. The (algebraic) dissi- 
pator A]j is a positive-definite matrix 1 . Note that the 
"dissipation" generated by the dissipator is that of states, 
including decoherence. In fact, this notion of "state dissi- 
pation" can be given a more precise geometrical meaning 
as follows. For any distance D on the space of density 
operators which is constructed from a monotonic metric 
(e.g. trace distance or Bures distance), completely posi- 
tive evolution cannot cause any Hilbert-space distances 
to expand [ !]: 

D[g(t) Pl ,g(t)p 2 ] <D[ Pl , P2 }. (8) 

From this result it is then easy to prove that positive- 
definite dissipators contract the state-space volume, 
whereas negative-definite dissipators expand the state- 
space volume (since they appear to be time-reversed con- 
tractions) . 

Derived in App. A and also with more context in 
Ref. [17], the second-order algebraic dissipator Ajj(t) 
evaluates to 

£ f dr [dr' (i|L m (r) a nm (r', t)(j\ L„(r') \j>) , 

nm J ° Ja 

(9) 

in the interaction (Dirac) picture. Expression (9) will 
be shown to be a positive-definite quadratic form for all 
microscopically derived noise correlations, thus agreeing 
with the Lindblad-GKS theorem in as much as is re- 
quired. Although the Markovian generators which ap- 
pear in the Lindblad equation are well known, to our 
knowledge these non-Markovian generators (which are 



The master equation for p exhibits the same structure as the 
right-hand side of Eq. (7). However, whereas in the Marko- 
vian regime one would simply have in place of & and A/j their 
time derivatives, this is no longer the case in the non-Markovian 
regime because <I> and # do not commute. 
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strictly algebraic and do not appear in the master equa- 
tion) are a novel discovery. 

The key to accomplishing our stated goal rests in the 
comparison of perturbative dissipators. We first note 
that from its microscopic origins, Eq. (2), the environ- 
ment correlation function is Hermitian in the sense of 

a{t,T) = a\r,t) , (10) 

where the boldface notation denotes a matrix with re- 
spect to the indices n, m in Eq. (2), and also positive 
definite in the sense of 

/ dr x f d7*ft(Ti)a(Ti,7a)f(7s) >0, (11) 
Jo Jo 

for all vector functions f(i) indexed by the environment 
correlator. All quantum correlations lead at least to 
accumulated decoherence since their algebraic Lindblad 
dissipator, Eq. (9), is necessarily positive definite. Ac- 
cumulated decoherence only implies that there is more 
net decoherent evolution than recoherent evolution. In 
general the stricter property of instantaneous decoher- 
ent evolution, A(i) > (boldface denotes here a matrix 
with indices /, J) , can only be satisfied by environements 
with local correlation function (Markovian processes) and 
would always produce a Lindblad master equation. How- 
ever, some very restricted classes of system-environment 
interactions, such as the interaction Hamiltonian in the 
rotating- wave approximation (RWA) [34], can be con- 
strained by the particular form of their coupling to be 
instantaneously decoherent. This characterizes the class 
of systems with non-Markovian dynamics (nonlocal en- 
vironment correlation function) whose master equation 
is, nevertheless, naturally of Lindblad form, though not 
necessarily at all times. 

The key result of this work is showing that the environ- 
ment correlation function itself provides a natural com- 
parison of state dissipation or decoherence strength. If 
two correlation functions are ordered c*+(i, r) > cc- (t, r) 
in the sense of positivity (11), then their correspond- 
ing second-order Lindblad dissipators are also ordered 
A+(t) > A_(£), and we can, therefore, say that one 
environment generates more state dissipation than the 
other, regardless of the system. For instance, the set 
of univariate Markov processes is totally ordered by the 
scalar magnitude of the respective delta correlations, e.g. 
2 S(t — t) > 1 S(t — t). In general, the set of all quantum 
correlations is only partially ordered, but nontrivial or- 
derings do exist. We illustrate this principle with several 
examples below. 

III. THERMAL CORRELATIONS 

A. Individual reservoirs 

Time-independent coupling to a thermal reservoir will 
always produce time-translation-invariant environment 



correlations which can be expressed in the Fourier do- 
main as 

a(w) = 7(w) [Mw) - w] . (12) 
%t(u)) = w coth(|^) , (13) 

in terms of the damping kernel 7(w) (anti-derivative of 
the dissipation kernel) and fluctuation-dissipation kernel 
kt(oj). This can be derived directly from first prin- 
ciples in Eq. (2), by demanding a coupling- invariant 
fluctuation-dissipation relation (FDR), or by demanding 
coupling-invariant detailed balance in the master equa- 
tion [17]. 

Note that a Markovian quantum regime (complex 
white noise) necessarily implies a local damping kernel 
and high temperature (local FDR kernel). The Marko- 
vian regime is reached when the system time scales are 
much slower than those of the environment, so that we 
can take the zero-frequency approximation 

Urn &(oj) = j(0) 2T . (14) 

u— >0 

Markovian processes can, therefore, be ordered in their 
decoherence strength by the product of their damping 
and temperature, a result which is well known. If one 
inquires as to the temperature of an unknown Markovian 
process, the FDR kernel (or noise-to-damping ratio) will 
always reveal this. 

In general, nonlocal correlations (e.g. for finite temper- 
ature), and thus decoherence strengths, are not totally 
ordered. For a fixed temperature thermal correlations 
can be ordered by damping. On the other hand we have 
the inequality 

KhotM > K co i d (w) > \oj\ , (15) 

and so for fixed damping, correlations can also be ordered 
by their temperature. Therefore, finite-temperature ther- 
mal correlations are partially ordered by damping and 
temperature. 

7strongH [KhotM - W] > Tweak H [Scold (w) ~ w] , 

(16) 

and so it immediately follows that 

d^ r t ong H>«weakM- (17) 

If one environment has weaker damping but a sufficiently 
higher temperature, then the two correlations cannot be 
ordered - the implication is that different systems would 
decohere faster or slower for each environment, but not 
in a manner which can be strictly ordered. 

From Eq. (16) we can now compare environments of 
low temperature and nonlocal damping. For fixed damp- 
ing, the monotonic ordering of temperature is no surprise. 
While for fixed temperature, the ordering of damping is 
more subtle though also not surprising. The damping can 
be increased by an overall rescaling, say 7(0; ) — > 2~y(uj), 
or multiplying it by a frequency-dependent function g(uj) 
such that g(uj) > 1, Vw. 
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B. Multiple environments 

Here we wish to recover and generalize the work of 
Beer & Lutz [7] wherein they compared the decoherence 
rates of collective environments with different tempera- 
tures and ohmic cutoff frequencies, specifically for linear 
coupling to an oscillator and with both reservoirs at a 
relatively high temperature. For multiple environments 
we can make the same comparison by using the natural 
measure of decoherence strength from the Lindblad dis- 
sipator and environment correlation function. First we 
note that for coupling to one reservoir, the individual 
thermal correlations can be expressed via Eq. (12). Next 
we note that for any monotonic cutoff regulator, with 
fixed local limit 7(0) and variable cutoff A, then 

7hi g hH > 7iowH , (18) 

where Ahi g h > Ai ow (referred to as "fast" and "slow" 
in Ref. [7]). We can also compare the individual FDR 
kernels as per Eq. (15). Finally we can use the above 
relations to construct the mathematical inequality 

[7hi g hH - 7iowH]{[^hotM -u]- [Scold (a>) - w]} > , 

(19) 

which can then be rearranged to show that 

o£&(o;) + > a£» + ag&M , (20) 

which is consistent with the results of Ref. [ ], when 
we interpret the left and right-hand sides of Eq. (20) 
as comparing the decoherence strengths of two different 
collective (non-equilibrium) environments. Note that as 
the individual reservoirs are Gaussian and independent, 
one may simply add their correlations in determining the 
collective correlation. Our result applies more generally 
than that of Ref. [7] , in terms of coupling and tempera- 
ture, though we do not calculate a specific decoherence 
time. Their work relied upon what is essentially the exact 
FDR kernel, but which has been referred to as an effec- 
tive temperature [35] in the classical regime. We would 
rather avoid this nomenclature given that such environ- 
ments will lead to an asymptotic stationary state which 
is not thermal in general (sufficiently simple systems can 
reach a thermal state, but the corresponding temperature 
will be different for different systems). 

IV. DISCUSSION 

In this work we have motivated a general notion of 
decoherence strength as that generated by the quantum 
correlations of the environment, which in turn determine 
the magnitude of the algebraic Lindblad dissipator and 
thus "state dissipation" itself (see Sec. II for its precise 
meaning) . The ordering of decoherence strengths in this 
formalism is only partial, although when it occurs it is 
independent of the initial state of the system and of the 
the particular operators L„ which characterize the sys- 
tem coupling to the environment. However, this is not to 



say that comparison of environment correlations is not 
useful when there is no strict ordering. 

State-dependent decoherence is of particular interest 
in the search for decoherence-free states and how they 
emerge in certain classes of environments. If for two en- 
vironments ol\(uj) — a.2{uj) is indefinite, then there could 
be states corresponding to the vectors f in Eq. (11) which 
exploit this. As an explicit example, the well-known phe- 
nomena of sub-radiant and super-radiant spontaneous 
emission of atoms in the electromagnetic field vacuum 
can be viewed in this manner. A detailed analysis can be 
found in Ref. [36], which is formulated most directly in 
terms of the electromagnetic field correlation function. 
A well-known result covered there is that for a pair of 
two-level atoms very close together there is a joint state 
which decoheres rapidly (super-radiance) and another 
joint state which decoheres very slowly (sub-radiance), 
as compared to the decoherence rate of an isolated atom. 
This is possible because when two atoms are brought 
close together in the electromagnetic field, their multi- 
variate correlation ei ncar (cj) cannot be totally ordered 
with respect to the factorized correlations present in the 
far- distance limit, Q!f ar (u;). 

As mentioned in the introduction and defined precisely 
in Sec. II, strictly speaking the ordering that we have 
established is for environment-induced irreversibility (in 
terms of "state dissipation"). For weak coupling, one 
expects that in most situations of interest decoherence 
will dominate over thermalization at short times in such 
an irreversibility process. In fact, if one considers times 
much shorter than the relaxation timescale(s) and the 
real part of the a(t, t) can be neglected, the decoherence 
strength for an environment in equilibrium is character- 
ized by 7(w)kt(w). Furthermore, from Eqs. (16) and 
(19) one can immediately see that the same inequalities 
as in Eqs. (17) and (20) also hold for ~f{u})Rr{w). 
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Appendix A: Non-Markovian dissipators 

As mentioned in Sec. II, application of the Lindblad- 
GKS theorem in the non-Markovian regime requires, not 
the time-translation generator, but the algebraic genera- 
tor. To obtain the non-Markovian dissipator quickly, we 
will first begin with the von Neumann equation for the 
density-matrix propagator of the closed combined system 
C, consisting of system + environment, in the interaction 
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(Dirac) picture: 



d „ / \ 



(Al) 



where Q c (t) is the interaction-picture propagator and 
C_i(t) is the interaction-picture Liouvillian, defined as 

g c (t) = g^(t)g c (t), (A2) 
c l (t) = g F 1 (t)c 1 (t)g F (t), (A3) 

in terms of the free (non-interacting) system + environ- 
ment propagator £7f(*)- Note that g_ G (t) are superop- 
erators, which act on usual Hilbert-space operators as 
g c (t)A = U(t, 0) A Ur x {t, 0), and similarly for Q F (t) 
and its inverse. The interaction Liouvillian is further- 
more defined in terms of the interaction Hamiltonian 
Hi(*) as 



C J (t)p = -i[H I (t),p] 



(A4) 



For simplicity we consider factorized initial states 
of the system and environment so that we can solve 
Eq. (Al) with a Neumann series and then trace out 
the environment more easily. (Properly correlated initial 
states of the system and environment within this formal- 
ism are considered in Ref. [37].) The Neumann series 
produces a perturbative expansion of the open-system 
propagator: 




ft -i oo 

/ drC Y (T) =1 + (A5) 
Jo -I k=i 



o 
n-i 



(A6) 



also in the interaction picture and where To = t. This 
series can be contracted into the single exponential 



g_(t) = , 
g(t) = g (t) e * (t) 



(A7) 
(A8) 



where Qo(t) is the free propagator for the system. For 
symmetric noise (with vanishing odd cumulants), such as 
Gaussian with vanishing mean, the perturbative genera- 
tors can then be found to be 



(A9) 
(A10) 



This is equivalent to solving the master equation via . 
nus series [38] in the interaction picture. It should be 



noted that Magnus-series solutions are slightly secular in 
time, since in general the Magnus series has a finite ra- 
dius of convergence [ ] . In this context the second-order 
Magnus-series solution will accurately match the correct 
solution to second order at early times, and then only 
match the correct solution to zeroth order at later times, 
wherein it converges to the RWA solution. For some as- 
pects of the solution this accuracy can be improved with 
a less-secular integrator. However, a careful analysis of 
the master equation and its solutions has shown that, 
due to unavoidable degeneracy, the second-order master 
equation is fundamentally incapable of providing the full 
second-order solutions [40] . Therefore we would not pro- 
mote these solutions as the second-order solutions, but 
they contain the most information pertaining to the non- 
Markovian dissipation, which one can extract from the 
second-order master equation. 

The Magnus-series solution to the second-order master 
equation gives rise to the second-order algebraic genera- 
tor 

*(*) = fdr fdr' (£ x (r) Cj(r')) E + 0{C\) , (All) 
Jo Jo 

in the interaction picture. In terms of the interaction 
Hamiltonian, the Lindblad coefficients of the algebraic 
generator 3?(i) are then 

A u (t) = fdr fdr'Um^r) (f\Ui(r') |i» E , 
Jo Jo 

(A12) 

given the representation in Eq. (7). With the interaction 
Hamiltonian expanded as a sum of separable operators, 
as per Eq. (1), the coefficients evaluate to 



dr dr' (i\ L m (r) a nm {r\ r) (j\ LJr') \j>) 



(A13) 

in terms of the environment correlation function. Both 
forms are positive-definite quadratic forms. This im- 
plies that the second-order master equation must gener- 
ate completely-positive maps to second order. Further- 
more, it also implies that the second-order Magnus-series 
solution g(t) = go(t) e^w is actually completely posi- 
tive exactly (rather than just through second order). 

Therefore, one can see that the correlation function 
a(£, r) is not only the influence kernel and complex noise 
kernel, but also a decoherence kernel which determines 
the magnitude of the non-Markovian dissipator, at least 
in the weak-coupling limit. Note that in this derivation 
the correlation function may describe multivariate noise 
arising from an environment completely out of equilib- 
rium. 
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